Light propagation in chiral media 
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Light propagation in chiral media is discussed. We derive the wave equations for a probe pulse 
propagating through a chiral medium, and solve them analytically in Fourier space using the slowly 
varying envelope approximation. Our analysis reveals the influence of the different medium response 
coefficients on the propagation dynamics. Applying these results to a specific example system, we 
show that chiral interactions already become important at experimentally accessible parameter 
ranges in dilute vapors. The chirality renders the propagation dynamics sensitive to the phase 
of the applied fields, and we show that this phase- dependence enables one to control the pulse 
evolution during its propagation through the medium. Our results demonstrate that the magnetic 
field component of a probe beam can crucially influence the system dynamics even if it couples to 
the medium only weakly. 

PACS numbers: 42.50.Gy,42.65.Sf,42.65.An 



I. INTRODUCTION 

Quantum optics describes the interaction of matter 
and light at the quantum level [U, In many cases, the 
underlying principles can already be understood from a 
fairly simple analysis, focusing on single atoms interact- 
ing with the electric field component of the driving laser 
fields. But recently, motivated by the aim of achievinga 
negative index of refraction in atomic media 
a qualitatively different class of setups received a con- 
siderable amount of interest. The key feature of these 
setups is that a probe beam interacts with the medium 
both with its electric and its magnetic field component 
at the same time. Only if both the magnetic and the 
electric response are suitable, a negative refractive index 
(NRI) can be achieved. Related to the need for a high 
magnetic response, these systems have two more distinct 
features in common. First, they operate at a high den- 
sity, such that the different atoms in the medium are 
no longer independent [1, [^, Second, in order to 

enable a certain enhancement mechanism for the mag- 
netic response, the laser fields couple to the medium in 
such a way that the system becomes chiral [ll|, [l^ . 
Unfortunately, an experimental implementation of these 
schemes to achieve atomic NRI media is very challeng- 
ing. For example, suitable atomic species are rare, and 
the required high density implies a number of problems 
such as Doppler broadening, radiation trapping [T3l.[l3|. 
dephasing or unwanted non-linear processes. 

It is important to note, however, that these challenges 
to a large degree arise from the specific requirements for 
NRI media. But the two key ingredients, high density 
and chirality, are of interest in their own rights. This 
raises the question whether it is possible to study the 
effects of density or chirality individually at much more 
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accessible parameter ranges. It turns out that the answer 
is affirmative. Regarding experiments at high density, 
potential model systems could be, e.g., Rydberg-atoms, 
doped crystals, or cold gases. In a previous work, we 
studied dipole-dipole interactions in an electromagnet- 
ically induced transparency setup at high density p^ . 
We could show that the polarization contribution to the 
local field strongly modulates the phase of a weak probe 
pulse. It was found that this phase modulation distinc- 
tively differs from the nonlinear self-phase modulation 
experienced by a strong pulse in a Kerr medium. 

Here, we turn to the second key ingredient of NRI 
schemes, and analyze light propagation in a dilute sam- 
ple of chiral atoms. Quantum optical setups routinely 
neglect the magnetic field components of the applied elec- 
tromagnetic fields altogether. Generally, this is justified 
due to the weak interaction strength, and since mag- 
netic transitions are typically off-resonant from the ap- 
plied fields. Interestingly, a complete switch from electric 
to magnetic dipole transitions would not change the ob- 
tained results qualitatively, as the structure of the system 
Hamiltonian is essentially unchanged. It is only the com- 
bination of magnetic and electric field interactions in a 
chiral setup that gives rise to different results. 

We start by solving the wave equations in slowly vary- 
ing envelope approximation (SVEA) for an arbitrary chi- 
ral medium. Our solution allows us to understand the 
influence of the different response coefficients that char- 
acterize a chiral medium on the light propagation. In a 
next step, we discuss the necessary conditions to imple- 
ment a chiral atomic medium. Based on these results, 
we analyze an example system with chiral structure. We 
study the light propagation through this medium, and 
find that chiral interactions become relevant already at 
low medium density. The chirality gives rise to a sensi- 
tivity of the system dynamics to the relative phase of the 
applied fields. We show that this dependence can be used 
to control the propagation dynamics of a probe pulse, 
even throughout the pulse propagation. The results are 
interpreted based on our analytical solution, and all con- 
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elusions are verified using a full numerical solution of the 
propagation equations. 

In a broader context, our results demonstrate that the 
magnetic field component of the probe field can crucially 
influence the propagation dynamics of the light pulse, 
even though it couples to the medium only very weakly. 
This may pave the way to further applications of mag- 
netic quantum optics, utilizing both electric and magnetic 
components of the applied laser fields. Also, it will turn 
out that the analysis of chiral media sheds light on the 
related so-called closed- loop systems [IB, [13 • We wfll 
find that the chiral media are an ideal implementation of 
a closed-loop phase-control scheme of light propagation. 

The paper is structured as follows. In Sec.|TTl we define 
and analyze chiral media in terms of parity transforma- 
tions, and discuss their linear response. The following 
Sec. mil analyzes light propagation in chiral media. In 
particular, we analytically solve the chiral wave equations 
in SVEA. The following Sec. IIVI discusses possible imple- 
mentations of chiral media in atomic systems. Based 
on these results, we discuss light propagation in a con- 
crete chiral model system in Sec.|Vl with the emphasis on 
phase control of the propagation dynamics. Section IVII 
discusses and summarizes the results. 



II. CHIRALITY 

An object is defined to be chiral if it is not invariant 
under parity transformation. Classical electrodynamics 
is known to be a parity invariant theory. To see this, and 
as a basis for our following analysis, we start from the two 
inhomogeneous macroscopic Maxwell's equations [18], 



V X E = -dtB , 

V X H = 9(0 , 



(la) 
(lb) 



with the electric field E, the magnetic field B, the elec- 
tric displacement field D, and the magnetizing field H. 
The medium is described by its polarization P and mag- 
netization M, which are related to the quantities entering 
the Maxwell's equations via the constitutive relations 



D = eoE + P, 



H 



1 
Mo 



-B-M. 



(2a) 
(2b) 



Due to the used SI unit system, the permittivity Sq 
and the permeability //q of free space occur. We do not 
consider the two homogeneous Maxwell's equations since 
with no external sources or currents and assuming trans- 
verse fields they are always fulfilled. 

Since E, P, and D are polar vectors they change sign 
under a parity transformation. In contrast, B, M and H 
are unchanged, because they are axial vectors. Together 
with the sign change of the curl operator, Maxwell's equa- 
tions Eqs. ([!]) as well as the constitutive relations Eqs. 
are parity invariant. 



While the governing equations are parity invariant, 
this does not need to be true for the macroscopic medium 
itself. Assuming a homogeneous isotropic medium, in 
linear response theory we can define linear electric and 
magnetic susceptibilities xe and xh, 

F{t)= J eoXE{T)-E{t~r)dT, (3a) 
M(t) = J Xff(r)H(t-r)dr. (3b) 

These equations only relate quantities of the same, odd 
or even parity and we see that the coefficients xe and 
Xh are scalars and thus unaffected by a parity trans- 
formation. From the perspective of an electrodynamic 
description, the corresponding medium is not chiral. A 
different situation arises, if a medium also exhibits cross 
couplings in a way that an electric field induces mag- 
netization and a magnetizing field induces polarization. 
In that case Eqs. ^ have to be extended by two extra 
terms, 



Pit) = J eoXEirMt - r) + ^^EHir)li{t - T)dT, 



(4a) 

M(i) = / XHiT)n{t - r) + -^^HE{T)E{t - T)dT, 

(4b) 

with the cross coupling coefficients ^eh, ^he, and c = 
1 / y^EQ^o the speed of light in vacuum. Under a parity 
transformation, the cross coupling coefficients (,eh and 
S,HE change sign, and thus are pseudo-scalars. If these 
cross couplings are non- vanishing, the medium properties 
change under parity transformation and we deal with a 
chiral medium. 

In Eqs. (gl), we have assumed that the medium and 
light field configuration are such that the response co- 
efficients are (pseudo-) scalar rather than tensorial as it 
is the case in general. It turns out that this is possible, 
e.g., if the probing light field is circularly polarized in 
the X ~ y plane and propagates in ±z direction [B, 0]. 
This configuration also allows for a different perspective 
on chirality. A circularly polarized light wave is a chiral 
object itself, since a left circularly polarized light wave 
under parity transformation is mapped onto a right cir- 
cularly polarized light wave and vice versa, rather than 
onto itself. Because of global parity invariance of electro- 
dynamics, considering a medium interacting with a right 
circularly polarized light wave is equivalent to consider- 
ing the inverted medium interacting with a left circularly 
polarized light wave. It follows that a medium is chiral if 
its response to a left and right circularly polarized light 
wave is different. This difference manifests itself in the 
refractive index, which in a chiral medium with response 
as in Eq. (jl]) was shown to evaluate to 



n± = xjefi- ^{^EH + ^he)^ T ^{^EH - ^he) ■ (5) 
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Here, e = 1 + xe and /z = 1 + xh- The upper sign applies 
for left circularly polarized light and the lower sign for 
right circularly polarized light. We see that in Eq. ([5]) it 
is equivalent to change the sign due to selecting a differ- 
ent polarization or to change the sign of ^eh and S^he 
due to inverting the medium. Nonzero cross coupling co- 
efficients £^EH and £^he in both cases reflect the chirality 
of the medium. 



III. LIGHT PROPAGATION IN CHIRAL MEDIA 

In this section, we derive the wave equations for the 
electric and magnetic field components of a probe field 
propagating through a chiral medium. This is in contrast 
to the standard treatment which assumes a coupling to 
the electric component only. We apply the slowly varying 
envelope approximation (SVEA) in order to transform 
the second order wave equations to equations involving 
first order derivatives only. Using Fourier transformation 
techniques, we find the general solution to the first order 
wave equations, and compare the result to the index of 
refraction obtained from a single particle susceptibility 
analysis. 



A. Wave equations for the electric and tiie 
magnetic field components 

From Eqs. ([T]) and we can derive wave equations 
for E and B with sources described by P and M as. 



equations for the envelope functions, 



A 



E = fiodfP + fiodt"^ X M, (6a) 
B = ^oAM - fiodtV X P. (6b) 



In the following, we specialize to the case of a one- 
dimensional propagation. We choose the positive z axis 
as the propagation direction and separate a complex en- 
velope function Xo{z, t) from the carrier wave. 



1 



X = -Xaiz,t)e±e 



i{kQZ—uJc)t) 



H- c.c. 



(7) 



where X stands for E, B, P, or M. The complex unit 
vector for left or right circular polarized light is e± , and 
the wave number of the carrier wave in vacuum is fco with 
ujo = cko as the corresponding frequency. Assuming that 
variations in space and time of the envelope function are 
on a much larger scale than the wavelength Aq = 27r/fco 
and the oscillation period To = In juj^ of the carrier wave, 
we can apply the slowly varying envelope approximation 
(SVEA) [1| . In this approximation, derivatives of the en- 
velope function in space and time are neglected compared 
to derivatives of the carrier wave. As a consequence, 
SVEA reduces the wave equations Eqs. (O to first order 



9. + -dt 



+ -dt 
c 



Eo{z,t)^'^Poiz,t)T,^Moiz,t), (8a) 
2£o 2£oc 

B,{z,t)^'^Moiz,t)±^P,{z,t). 

(8b) 



As before and in the following, the upper sign applies for 
left circular polarization whereas the lower sign applies 
for right circular polarization. Eqs. ([5]) are considerably 
simpler to solve analytically as well as numerically than 
the second order Eqs. ([6]) we derived initially. Yet, they 
still incorporate the essential physics of light propagation 
in one spatial dimension for a chiral medium with a cou- 
pling to both the electric and the magnetic component 
of a weak probe field. 



B. Solution of the wave equations 

We will now solve the wave equations ^ in SVEA 
to a level that enables us to determine how the actual 
propagation dynamics depends on the four quantities xe, 
Xh, ^eh, and ^he characterizing the medium. First, 
we transform Eqs. ^ into Fourier space. To substitute 
Pq and Mq, we also Fourier transform Eqs. Next, 
to replace Hq by Bq, we use the Fourier transform of 
Eq. (j2bp . Finally, we can combine both wave equations 
into a single matrix equation for the vector F(z, Ap) = 
{Eo{z,Ap),cBoiz,Ap))^, 



9,F(z,Ap) =ifcoA^F(0,Ap), 
and the elements of the matrix M are given by 



(9) 



Mi^l^ ^ + ^{XE~^EHiHE±i^HE), (10a) 



Ml,2 = ■^^{S,EH ±iXH) , 

M2,l = ^^i^HE T iMXB ± iS.EHS.HE) 

M2.2 = ^ + ^{XH ± i^EH) ■ 

UJQ Zfl 



(10b) 
(10c) 
(lOd) 



In these equations, the detuning Ap — uj — ujq accounts 
for the frequency distribution of the envelope functions 
around the carrier frequency. Applying the vacuum 
phase relation for circularly polarized light, cBq — ^iEq, 
we find that the initial condition for the evolution Eq. 



F(OA)-(' 



(11) 



is an eigenvector to matrix M. in Eq. Therefore, 
the solution can be calculated with the help of the cor- 
responding eigenvalue 77. The solution reads 



F(z,Ap) 



£;o(0,Ap) 
TiEo{0,Ap) 



(12) 
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and the frequency dependence is contained in the eigen- 
value 



[XE + - ^EH^HE T i{^EH - ^He)] 



(13) 



This resuh can be compared to the refractive index 
Eq. ([5|) in a chiral medium obtained from a single- 
particle susceptibility analysis. The application of SVEA 
is equivalent to the condition \n± — 1| <^ 1, such that wc 
can perform a first order Taylor expansion of the square 
root contribution in Eq. ^ . Keeping terms linear in the 
response coefficients we find 



1 \ * 

n± = l + -[XE + Xh) T i^{.£.eh 



i 



HE) 



(14) 



Linearizing Eq. (|13p and taking into account the contri- 
bution of the carrier wave, we find the same result as 
in Eq. (fT4|). This demonstrates the consistency of our 
solution with the previously calculated chiral index of 
refraction. Our result based on the solution of the wave 
equations, however, has the advantage that it offers the 
possibility to study propagation dynamics. Because of 
the simple formulas we have derived, analytical solutions 
become accessible. If for a specific system the frequency 
dependencies of xe, Xh, £,eh, and £_he are known, they 
can be used together with the solution Eq. p2)) in fre- 
quency space to calculate the actual dynamics in the time 
domain. 



IV. CHIRAL ATOMIC MEDIA 

In this section we discuss the general conditions for 
realizing a chiral atomic medium in the form as intro- 
duced in Sec. [Ill As a first condition, a medium in addi- 
tion to an electric probe field transition should exhibit a 
near-degenerate transition coupling to the magnetic com- 
ponent of the probe field. A second condition is that, to 
enable chiral cross couplings, there must be a mechanism 
to allow for polarization [magnetization] to be induced by 
the magnetic [electric] probe field component. 

The first condition can be difficult to fulfill in atomic 
media. A small frequency difference between electric and 
magnetic probe field transition could be removed by an 
external magnetic field leading to a Zeeman shift, but 
magnetic dipole transitions are usually found in a much 
smaller frequency range compared to electric dipole tran- 
sitions. The reason is that due to selection rules imposed 
by the corresponding transition matrix element, mag- 
netic dipole transitions can only occur between states 
of the same parity. This implies that magnetic dipole 
transitions only take place between states with the same 
angular momentum quantum number whereas electric 
dipole transitions connect states of different angular mo- 
mentum. A possible solution are atomic species with a 
large spin-orbit coupling [I^ . Then, states with the same 



angular momentum quantum number can display a large 
splitting and their transition frequency can become com- 
parable to the transition frequency of an electric dipole 
transition. This could also help to realize a probe field 
transition frequency in the optical domain which is es- 
pecially attractive for applications. Motivated by ideas 
like NRI, several atomic species have been investigated 
with respect to this condition. Dysprosium Hy- 
drogen 4] , metastable Neon and other noble gases J, l2l| 
have been identified as promising candidates. An alter- 
native approach could be to make use of non-resonant 
two-photon transitions in order to relax the requirement 
of having both electric and rnagnetic transitions at simi- 
lar transition frequencies 0, l2l| . 

The second condition can be implemented in closed- 
loop systems. In these systems, the laser-driven transi- 
tions form a closed interaction contour such that pho- 
ton emission and absorption can take place in a cy- 
cle [3 1221 . Because the closed interaction contour en- 
ables quantum mechanical pathway interference, it can 
render the system phase dependent. Closed-loop sys- 
tems have mainly been investigated in non-chiral con- 
texts [l3,[2l,l2ll2i,[2i,|2l[28]. Here, the closed interac- 
tion contour only includes the coupling to one (electric) 
probe field component. This is known to induce extra 
contributions to the medium response, but these cannot 
be utilized for light propagation of a pulsed probe field 
in a straightforward way [23, [13] • The reason is that al- 
ready without the probe field, such a configuration would 
give rise to parametric processes which scatter the control 
fields into the probe field mode [H, [30|, IM HI, HI . As 
we will see in Sec. |Vl this conclusion changes in the case 
of chiral closed-loop systems. There, we will find that 
the closed-loop phase can provide a convenient control 
parameter for light propagation in chiral media. 



V. CHIRAL CONTROL OF PROPAGATION 
DYNAMICS 

A. Model system 

In this Section, we apply the results obtained so far 
to discuss light propagation in a specific system with 
simultaneous coupling to the electric and the magnetic 
component of a probe laser field. The atomic medium is 
modeled as the five level system shown in Fig. [1] It was 
first introduced in the context of generatingatomic me- 
dia exhibiting a negative index of refraction [6] . In zeroth 
order of the probe field, the atomic medium forms an ef- 
fective three-level system in A-configuration consisting of 
states ]1), ]4) and ]5). Two strong resonant control fields 
with Rabi frequencies fJi and ^2 prepare the atoms in a 
dark state, i.e., a stable coherent superposition of states 
]1) and ]4). The magnetic and electric probe field compo- 
nents couple to degenerate magnetic and electric dipole 
transitions with Rabi frequencies and I^e, and the up- 
per levels of these two dipole transitions are coupled by 
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an additional resonant control field with Rabi frequency 
ilc- We assume that state |1) and |4) have the same 
parity. Because the electric probe field transition cou- 
ples states of different parity whereas the magnetic probe 
field transition couples states of the same parity, the addi- 
tional control field couples to an electric dipole transition. 
Note that the Rabi frequencies ili — \ exp[i{kif+ ipi)] 
are complex, and contain the wave vector ki and the ab- 
solute phase ipi {i £ {1, 2, C, E, B}). 

As discussed in Sec. HI the linear response of the chi- 
ral atomic medium is characterized by two contributions. 
First, direct contributions arise from the polarization 
[magnetization] induced by the electric [magnetic] probe 
field component. It turns out that the direct contribution 
to the magnetic response is typically small. The reason is 
that magnetic dipole couplings usually are much weaker 
than electric dipole couplings. 



|3). 



m/d ' 



ac 



(15) 



with magnetic dipole moment m, electric dipole moment 
d, and fine structure constant a. Therefore, in compar- 
ison to the direct electric response, the direct magnetic 
response is parametrically suppressed by two powers of 
the fine structure constant 1341. 



XH 



(16) 



The other contributions to the medium response are cross 
couplings, which arise due to polarization [magnetization] 
induced by the magnetic [electric] probe field component. 
It was found that the considered level scheme displays a 
relevant magnetic response due to these cross contribu- 
tions [s^ ]. A magnetic response arising from such a 
cross coupling is only suppressed by a factor on the or- 
der of a. Furthermore, it is proportional to the strength 
of the ground state coherence of state |1) and |4) and the 
control field Vic- Since cross couplings of this type lead 
to nonzero values for S^eh and S,hEi they also render the 
medium chiral. 

The Hamiltonian of the system in Fig. [1] in a suitable 
interaction picture is given by [l| 



h 

~ 2 

f^i A51 -I- 0.2 + H-c. 



■ {q.b A21 + Oe A34 + Oc A^2 



(17) 



where Ap is the probe field detuning and atomic transi- 
tion operator is defined as Ajk = An interesting 
insight into the involved physics can be gained by apply- 
ing a further unitary transformation to the Hamiltonian 



Hi = -h\{A22+ A^:i) 
_ h 
~ 2 

+ \ni\A^i + [f^al A54 



(|f}B|A2i + |fi£;|A34 + |r!c|e'*^32 



H.c. 



(18) 



n V 



fir 



|5)- 



" O. / / 



n 



X If- 



FIG. 1: (Color online) Atomic level scheme used to describe 
the chiral medium. Two strong control fields with Rabi fre- 
quencies rJi and VL2 prepare a coherent superposition of states 
jl) and |4). The magnetic and electric probe beam compo- 
nents couple to two transitions with Rabi frequencies Q, b and 
Q.E, respectively, fie describes an additional control field. 
The relative state parities are denoted by -I- or — . 



It can be seen that Hi contains a phase contribution 



= (}P2 - ^\ + ^c) + (h - fci + kc)r, 



(19) 



which cannot be removed entirely via unitary transforma- 
tions. The origin of this phase dependence is the closed- 
loop structure of the considered level scheme. Starting 
from |1), the system can evolve in a non-trivial loop se- 
quence |1) ^ |5) ^ |4) ^ |3) ^ |2) ^ |1) back to 
the initial state. This enables pathway interference, e.g., 
from |1) to |3) either via |5) and |4), or via |2). The 
phase difference between these two interfering pathways 
is equivalent to the phase difference given by Eqs. fTB]) 
and p^ . It is interesting to note that the closed- loop 
phase $ does not contain properties of the probe field, 
which is due to the fact that a closed transition loop in- 
volves both an absorption and an emission of a probe 
field photon. Furthermore, $ does not depend on time, 
because of the cancellation of the probe field and be- 
cause the three control fields ni,^l2,^c are applied on 
resonance. 

We are now also in the position to understand why 
a phase-control of light propagation is possible in chiral 
closed-loop media, in contrast to closed-loop media cou- 
pling only to the electric component of the probe field. 
From Fig. [l] it can be seen that in the presence of the 
probe field, a closed loop is established, such that the 
medium becomes sensitive to the closed- loop phase <I>. 
But in the absence of the probe field, two transitions 
in the closed loop structure are undriven. It follows that 
parametric processes scattering the control fields into the 
probe field mode cannot occur, independent of the phase 
matching condition assumed for the wave vectors of the 
probe and driving fields. This crucial difference to closed- 
loop media coupling to one (electric) probe field compo- 
nent only enables the phase control of light propagation 
that we will find in Sec. IV Bl 
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In the following, we continue our analysis based on the 
Hamiltonian Hi in Eq. pT]) and assume that the laser 
configuration satisfies the usual phase matching condi- 
tion k2 — ki + kc , such that the system dynamics only 
depends on the relative phase 



(20) 



In zeroth order of the probe fields, the equations of mo- 
tion (EOM) governing the dynamics of the atomic de- 
grees of freedom can easily be solved. The non- vanishing 
density matrix elements are 



(0) 

P44 



(0) 



(21a) 
(21b) 
(21c) 



The linear response to a weak probe field does not disturb 
this dark state prepared by fli and fl2, and can thus be 
calculated on the basis of the lowest order populations. 
In order to make connection to NRI calculations, we also 
take into account local field corrections (LFC) as they 
occur in denser gases, even though we will find later that 
magnetic field effects can already be observed at densi- 
ties low enough to neglect LFC. In the literature, a gen- 
eralized form of the Clausius-Mossotti equation [18[ has 
been used to numerically extract xe, Xh, ^eh, and £,he 
from the steady state solution to linear order in the local 
fields E'loc and Bioc in Hi • We pursue a different strategy 
here which provides us with an analytic solution includ- 
ing LFC [la|. We directly replace the local fields in the 
EOM by the external fields i?ext and i?oxt with the help 
of the Lorenz-Lorentz relations UM 



E\oc — Ecxt + 7^ P , 

3eo 



-Bloc — Be 



Mo 



M. 



(22a) 
(22b) 



This leads to nonlinear EOM if the polarization P and 
the magnetization M are expressed by microscopic expec- 
tation values of the electric and magnetic dipole density 
containing elements of the atomic density matrix. Since 
we are interested in linear response, we expand the EOM 
up to linear order in the probe field, and solve for the 
steady state. From this steady state solution, we find for 
the response coefficients 



Xe 



Xh 



^EH 



24 



t 34^ 24 H 4 — 



(23a) 



r2ir 



31 



-(■ 



iEH\p^!h 



3api?fr, 



31 



3L734Q; \p^^\ 



\£tc\ 
4 



r34r24 + 



(23b) 
(23c) 



£,HE 



3L734a 



(o)| Hd 

4 



r2ir 



31 



(23d) 



where L = /Att^ includes the dependence on the 
number of particles NX'^ per wavelength A cubed. The 
decay rate and the detuning including the typical LFC 
shift for the electric and magnetic probe field transitions 
are given by 



r34 - 7/2 + 7dcc - «(Ap + ^14^^734 /2) 

r2i = 721/2 + 7dcc - i\ , 



(24a) 
(24b) 



where 7 = 731 -f 732 +734 is the overall decay rate of state 
1 3), the linewidth of the magnetic probe field transition is 
721, and 7dec accounts for additional decoherence induced 
by atomic collisions. 

Interpreting the dark state of |1) and |4) as a single 
quantum state characterized by the coherence between 
the two states, both probe field transitions can be viewed 
as being part of a three level electromagnetically induced 
transparency (EIT) system [H, . The electric probe 
field and the control field Vic form a A-system for which 



r24 = 721/2 - lAy 



(25) 



contains the ground state decoherence and the two- 
photon detuning. No additional decoherence 7doc be- 
sides the decay of the magnetic probe field transition has 
to be included if we assume that levels |2) and |4) are 
nearly degenerate, as then phase-perturbing collisions in 
the medium lead to random, but correlated phase shifts 
of the two states @ . The magnetic probe field and the 
control field Vic form a ladder system in which the de- 
coherence and two-photon detuning relevant for EIT is 
described by 



Fsi = 7/2 + 27dcc - «Ay 



(26) 



By inspection of the structure of the solution for xe 
and xh^ we find that the direct electric and magnetic re- 
sponses of the EIT system are amended by LFC. The 
electric response features the usual EIT transparency 
window which is only degraded by the weak ground state 
decoherence induced by 721. The transparency window 
of the direct magnetic response on the other hand is not 
as prominent due to strong decoherence 7/2 -f 27dcc in 
the ladder system. But this does not affect the propaga- 
tion dynamics much since xh is suppressed by a factor 
compared to xe as discussed above. In particular 
for propagation dynamics at densities corresponding to 
L < 1, it can be neglected. On the contrary, the chi- 
ral cross couplings are only suppressed by a factor of a 
and do influence the propagation dynamics. As expected, 
they are proportional to \pfi\^ \Vlc\-, and depend on the 
relative closed-loop phase $0. At first glance it might 
look like a phase dependence is also induced to xh by 
LFC because it depends on $0. But a closer look re- 
veals that this phase dependence is exactly canceled by 
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FIG. 2: (color online) Propagation dynamics in a chiral 
medium. The closed-loop phase $o = tt/S is kept constant, 
see inset. The absolute value (solid black lines), the real 
part (dashed blue lines), and the imaginary part (dotted red 
lines) of the probe pulse envelope are plotted for successive 
propagation times. Arbitrary, but consistent units are used 
and successive plots are vertically offset by the overall phase 
at the center of the pulse. The spatial phase dependence 
of each pulse (dashed green straight lines) is compared to 
the value expected from the analytical solution (thick gray 
straight line). 



the factor (,eh- From these solutions we expect a prop- 
agation with a reduced group velocity and with a phase 
dependent refractive index induced by the chiral cross 
couplings. 

One further observation about the chirahty of the 
medium can be made. In the sense of a macroscopic elec- 
trodynamic description, the control fields and ^2 

are part of the medium. Thus, they change sign when the 
medium is inverted and lead to the expected chiral be- 
havior of i,EH and (,H_E as described in Sec. [ill 

With the analytic solution for XHi S,eh, and ^he 
at hand, we can now proceed to study propagation ef- 
fects in a chiral medium, focusing on control of the prop- 
agation dynamics via the closed-loop phase as concrete 
application. 



B. Phase switching 

In this Section we study the propagation dynamics of 
a pulsed probe field in the chiral medium discussed in 
Sec. IV Al Due to the direct electric response, the probe 
pulse propagates with a reduced group velocity, as it is 
the case in a usual EIT system. If we write the group 
velocity 



the group index Ug it is given by 



(28) 



1 + n„ 



(27) 



This is the standard expression in a EIT system, ex- 
cept for the scaling by the population fraction pf^ in the 
ground state of the electric probe field transition. The 
direct magnetic response can be neglected and all addi- 
tional influence of the medium on a probe pulse is due 
to the chiral cross coupling of the electric and magnetic 
probe field components. 

The chiral cross coupling is mediated by the closed 
interaction loop and depends on the relative phase $0 
between the three control fields. Since the probe field 
propagates with strongly reduced group velocity whereas 
the control field flc propagates at speed of light, a phase 
change of <I>o during the propagation of the probe pulse 
changes the frequency dependent refractive index for the 
probe field nearly instantaneously. In a similar fashion as 
a light pulse can be slowed down, stored, and retrieved 
by changing the amplitude of the control field [13], the 
phase of the control field can be used as a tool to control 
the refractive index for the probe pulse. 

To clarify this idea further, in the following, we present 
two exemplary results from a numerical solution of the 
atomic EOM coupled to the wave equations for the two 
probe field components and the control field. We assume 
a medium with a density corresponding to L — 0.01. For 
A — 795nm, this translates to iV ~ 8 x lO^^cm""^, which 
is low enough to neglect LEG. The control field strength 
is chosen as \^c \ — 27, and the dark state as symmetric, 
1^1 1 = \^2\, such that pf-^ = pf^ — 0.5 and ^4"^ = 
—0.5. In both cases, we propagate a Gaussian probe 
pulse with a width in time of cr = 5O/7. Its spectrum 
is well contained in the EIT transparency window such 
that it suffers very little absorption and broadening. We 
use a coUisional decoherence of jdec — O.57 and for the 
probe field we assume left circular polarization. 

In the first example, the phase $0 is chosen constant as 
7t/2. In Fig. [21 we show the corresponding propagation 
dynamics. The envelope of the electric probe field com- 
ponent is depicted for four successive propagation times, 
together with the real and imaginary part of the enve- 
lope and the pulse phase. The group velocity observed 
in the numerical calculation is consistent with Eqs. (|77|) 
and ([28)1 . We find that throughout the propagation, the 
phase of the pulse increases linearly. Thus, the real and 
imaginary parts of the envelope interchange while the 
absolute value of the envelope remains constant. This 
dynamics is due to a positive real index of refraction at 
resonance. Its theoretical value can be deduced from the 
solutions given in Sees. IIIIB] and [Vl Neglecting xh, the 
eigenvalue Eq. (|13p reduces to 

V=— + UxET^{iEH-iHE)]. (29) 

luq 2 

The term linear in Ap arising from xe leads to the re- 
duced group velocity given in Eq. (|27p . We approximate 
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FIG. 3: (color online) Propagation dynamics in a chiral 
medium with phase control. The closed-loop phase is switched 
from 7r/2 to zero and back to 7r/2 throughout the propagation, 
as shown in the inset. The absolute value (black solid lines) 
of the probe pulse envelope is plotted for successive propa- 
gation times. The scaling and offsets are used as in Fig. [21 
At the beginning and at the end of the propagation the spa- 
tial phase dependence of each pulse (dashed green straight 
lines) is compared to the value expected from the analytical 
solution (thick gray line) . During the time period of zero con- 
trol phase, probe pulse envelope phase does not change and 
the pulse experiences gain instead, as indicated by the thick 
dashed red line. 



S^EH and S^HE by the constant term at resonance lead- 
ing to a constant refractive index. Fourier transforming 
Eq. (jl2[) back into time domain, we find 

^e{z, t)\ _ (VLE{0,t- Z/Vg)\ ,[3koZ (nn.^ 

nB{z,t)J - [nB{o,t^z/vg)J ^ ' ^^'^^ 

in which the polarization-dependent (3 is given by 



27dec7 + 87L+l^c|= 



(31) 



where we neglected the radiative decay of the magnetic 
probe field transition 721 which is small compared to jdec- 

Inserting the parameters chosen in Fig. [21 we obtain 
a positive real value of (3 such that the phase of the en- 
velope increases constantly with propagation distance. 
In Fig. [21 we compare this phase dependence obtained 
from our approximate analytical model with the phase 
from the numerical propagation of the pulse at different 
propagation distances. At the center of the pulse the 
agreement between numerical and analytical solution is 
very good. Away from the center there are small devi- 
ations that increase with propagation distance. These 
deviations can be attributed to the neglected frequency 
dependence of S,eh and ^he- 

In a second example, the phase $0 is changed during 
the propagation. Initially, it is set to tt/2 as in the first 



example. During the course of propagation, it is switched 
to zero for an intermediate period of time and finally, it 
is switched back to tt/2. The numerical results for the 
corresponding propagation dynamics are shown in Fig. [3l 
and the time-dependent value of the phase is shown as 
inset in Fig. [3l In the initial and final part of the prop- 
agation with phase 7r/2, we find the same propagation 
dynamics as in Fig. [21 with reduced group velocity Vg 
and with a linear increase of the phase. But during the 
period with control phase switched to zero, the phase of 
the pulse envelope remains approximately constant. In- 
stead, its amplitude increases exponentially. As in the 
first example, this dynamics can be explained based on 
the index of refraction at resonance. Calculating f3 for 
$0 = in the analytical solution, we find a negative 
imaginary value corresponding to gain. Again, we com- 
pare the analytical solution to the numerical data. The 
phase dependence is given by a linear increase with a 
plateau in between, whereas the increase in amplitude 
follows an exponential function, see Fig. [31 Both analyt- 
ical solutions agree reasonably well with the numerical 
data, given that the finite switch period with the corre- 
sponding transient dynamics is not taken into account in 
the analytical calculation. 



These examples show that by changing the relative 
control field phase $0, the dynamics of the probe field can 
be substantially infiuenced throughout the propagation. 
In the second example, the refractive index experienced 
by a pulsed probe field with left circular polarization was 
changed from a positive real value n > 1 to an imaginary 
value representing gain and back. With different values 
for the control phase, also n < 1, absorption, or inter- 
mediate cases are possible. This is illustrated in Fig. IH 
which depicts the real and the imaginary part of f3 ob- 
tained from our analytical calculation against the control 
phase. 
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VI. DISCUSSION AND SUMMARY 

In summary, we discussed light propagation dynamics 
in a chiral atomic medium. The chirahty arises from a 
cross couphng of the electric and the magnetic compo- 
nent of a probe field induced by the medium. First, we 
derived wave equations for both probe field components 
in the chiral medium. Using the slowly varying enve- 
lope approximation, wc then solved these equations in 
Fourier space. Our solution elucidates the actual depen- 
dence of the propagation dynamics on the direct and chi- 
ral response functions in a macroscopic electrodynamic 
description. Based on these results, we discussed the 
general conditions necessary for realizing a chiral atomic 
medium. In a concrete example, we analytically deter- 
mined the response functions for a specific chiral model 
system from the microscopic quantum mechanical equa- 
tions of motion. As our main result, we demonstrated 



that chiral effects can crucially influence the propagation 
dynamics already at experimentally accessible parameter 
ranges. In particular, we find that chiral effects can al- 
ready be observed in dilute vapors. As an application, 
we showed that the chiral couplings render the medium 
sensitive to the relative phase of the applied fields. By 
changing the relative phase of one of the applied cou- 
pling fields, the propagation of a probe pulse can be 
controlled dynamically during its passage through the 
medium. From this, we concluded that chiral media are 
an ideal implementation of closed-loop phase control of 
light propagation. In a broader context, it is remark- 
able that the magnetic field component of the probe field 
can substantially modify the probe pulse propagation dy- 
namics, despite the weak coupling of the magnetic field to 
the medium. It remains to be seen if such magnetic field 
couplings may also find other applications in quantum 
optics. 
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